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Abstract 


The Maloney [1] and Luebbers [2] subcell methods can be used to implement dielectric or conducting 
sheets thinner than the cell size in the FDTD grid. The two methods differ in the treatment of the field 
components situated half a cell from the thin sheet. For a sheet of infinite size stricken by an incident 
plane wave, it has been shown by theory and experiments that the Luebbers method suffers from 
discrepancy at the wide incidence angle, while the Maloney method is exact at any incidence angle [3]. 
The comparison of Maloney and Luebbers methods has been later extended to dispersive sheets, more 
specifically, sheets made with the Debye media used in Bio-electromagnetism to model the Human 
Body tissues (A. Alkandari et al, AP-S/USNC-URSI Symposium, FR-UB.5P, Boston, 2018). 


This paper presents the further comparison of the two aforesaid methods in more realistic situations 
with objects of finite size composed of several thin sheets such as dielectric, conducting or Debye 
medium sheets. Comparisons with objects composed with both PEC surfaces and thin sheets are also 
presented. The results from numerical experiments show that the difference between the Luebbers 
and Maloney methods is affected by the problem to be solved, especially the geometry and size of the 
thin sheets, the physical parameters of the sheets (permittivity, conductivity, relaxation time), and such 
other parameters as the location of the observation, for instance, inside or outside an enclosure. The 
interest of using the Maloney method, which is accurate in any situation, but is more complex to 
implement in a general purpose computer code, is discussed, with a special emphasize on the case of 
sheets made with Human Body media. 


Index Terms. FDTD, Subcell method, Debye media, Human Body. 





Introduction 


a Layers of physical thickness << skin depth and wavelength. 


L1] Luebbers method and Maloney method. Maloney method is 
an addition to Luebbers method. It introduces a special treatment 
of the E field perpendicular to the layer. 


L Which method is the best? Is Maloney significantly better in 
bioelectromagnetics. Is it really needed? Its implementation more 
complex in Bio-FDTD codes. 


A] In this presentation 


. Extension of Maloney method to dispersive Debye media (Human 
Body media). 

. Comparison of methods with infinite in size layers and with realistic 
layers of finite size. 








The thin layer (or sheet) methods 


With both Luebbers and Maloney, £F 
parallel within the layer 1s advanced using 
the regular FDTD equation with fictituous 


In addition, Maloney split E perpendicular, situated ^ cell outside 
the layer, into two fields E ., and E ,,,, advanced as in vacuum and as in 
the layer (&, , o.), respectively. A combination of these two fields is then 
used in the advance of the H field in place of E, by replacing EF with a 


E ,, + BE ou where B = d/Ax and a = 1-P. 


out 
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Agreement of Maloney with analytical while Luebbers 1s erroneous. 


-» The Maloney method can predict exact value of Brewster angle 


tan(0) 2 Je, 
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Conducting Sheet 


Sheet 1 mm in thickness in 5-cm FDTD grid. Conductivity o,= 0.1 S (skin depth = 5 cm) 
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Luebbers erroneous at high frequency and at wide incidence angle 
while Maloney agrees with analytic reflection at any incidence and 
frequency. FDTD Experiment agrees with FDTD theory. 


7/13 





Extension of Maloney to Debye media 


Ex out Ex. a 


2D TM case Hz in 
Additional Ex in 


and Hz in _ 
e Hz outl ag eee With Maloney EXE 
COREOUPTOE = O s mue split on one side of the. 
advance of I layer. Valid only if media 1 
Hz out1 Faraday and 2 are identical. 
Law contour for O 
advance of 


Hz out We split E orma, on both 


S sides 
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Thickness d 
Debye medium 1 Debye medium 2 


1-For advance of Ex in we use Hz in » (Hz out14Hz out2)/2 (valid because Hz continuous at interfaces) 


2-Hz_out1 and Hz out2 advanced by Faraday contours that hold Ex in and Ex out 
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2. AEy AEx out AEx in where 
Ay outi = H; outi + Ax TAU aI = 
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- Agreement of Maloney, including a « pseudo-Brewster Angle » 


- Luebbers layer is erroneous at wide angles 


Debye sheet between two Debye media 


Sheet between two 
different Debye media 
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Incidence Anele (?) Incidence Angle (°) 


Luebbers and Maloney yield almost identical results 


Debye sheet between two Debye media 


Same observation with dielectric media: even with small ratio of permittivity, 
Luebbers close to Maloney (and analytical). 
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In Human Body, thin sheets are always situated between different media 
> Is Maloney needed in Bioelectromagnetics - Human Body - calculations ? 


-? An exhaustive study is currently performed to answer this question (we think 
the answer is Yes) 
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Finite size layer (one size is finite => 2D) 
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Reference calculation Ax= 0.2 mm 
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Thin layer calculations Ax = 1 mm 
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A large sheet of finite size (3D problem) 


incident wave 
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- During clear time where effect of boundary of layer not arrived: significant 


difference between Luebbers and Maloney. 


- After clear time: no difference between Luebbers and Maloney -» scattering by 


boundary preponderant. 


A small sheet (3D problem) 
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A dielectric box 
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— No difference between Luebbers and Maloney in the interior 
— Small difference outside the box. 


A conducting box 
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Again, Luebbers and Maloney superimposed in the interior, and very small 
difference outside the box (only visile with lower conductivity when 
reflection is very small). 


Conclusion 


L] Maloney technique always yields correct reflection and transmission of plane 
waves for dielectric sheets, conducting sheets, Debye sheets. 


L] Luebbers technique is simpler but reflection/tranmission of plane waves is 
erroneous at wide incidence angles. 


L] However, with objects (1.e. layers or sheets) of finite and small sizes, the 
Luebbers error seems negligible because the preponderant energy 1s the energy 
scattered by the edges. 


1L] If sheets form a cavity, Luebbers and Maloney fields in the interior are very 
close to each other. 


L1] In Human Body, it seems (to be confrmed) that Maloney is not needed because 
media on the two sides of the sheets are different. 


_] Finally, the need of implementing the Maloney depends on the problem. It 
seems in many case, especially in Human Body, Maloney is not needed. This is 
rather good news since Luebbers 1s significantly simpler to implement. sig 





References 


[1] G. Maloney and G.S. Smith, “The efficient modeling of thin material sheets in 
the finite-difference time-domain (FDTD) method”, /EEE Trans. Ant. Propag., vol. 
40, pp. 323-330, 1992. 


[2] R.J. Luebbers and K. Kunz, “FDTD modeling of thin impedance sheets”, IEEE 
Trans. Ant. Propag. Vol. 40, pp. 349-351, 1992. 


[3] J.-P. Berenger and F. Costen, “Theoretical comparison of methods for 
modeling thin dielectric or conducting sheets in the FDTD Grid”, IEEE Trans. Ant. 
Propag., vol. 67, Jan. 2019. 





